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We consider generalizations of the Kardar-Parisi-Zhang equation that accomodate spatial 
anisotropies and the coupled evolution of several fields, and focus on their symmetries and non- 
perturbative properties. In particular, we derive generalized fluctuation-dissipation conditions on 
the form of the (non-linear) equations for the realization of a Gaussian probability density of the 
fields in the steady state. For the amorphous growth of a single height field in one dimension we give 
a general class of equations with exactly calculable (Gaussian and more complicated) steady states. 
In two dimensions, we show that any anisotropic system evolves on long time and length scales 
either to the usual isotropic strong coupling regime or to a linear-like fixed point associated with a 
hidden symmetry. Similar results are derived for textural growth equations that couple the height 
field with additional order parameters which fluctuate on the growing surface. In this context, we 
propose phenomenological equations for the growth of a crystalline material, where the height field 
interacts with lattice distortions, and identify two special cases that obtain Gaussian steady states. 
In the first case compression modes influence growth and are advected by height fluctuations, while 
in the second case it is the density of dislocations that couples with the height. 



I. INTRODUCTION 

Non-linear stochastic partial differential equations ap- 
pear extensively in problems of equilibrium and non- 
equilibrium statistical physics. For systems in thermal 
equilibrium, the form of these equations is constrained by 
fluctuation-dissipation conditions 1] that ensure conver- 
gence of the steady state probability distribution to the 
appropriate Boltzmann weight . Non-equilibrium sys- 
tems are not similarly constrained, and there is no sim- 
ple way of finding their behavior in steady state (if any). 
However, there are examples in which steady states can 
be found exactly as solutions of the associated Fokker- 
Planck equations. In this paper we review some such 
examples, and introduce several new ones, along the way 
seeking general principles for finding steady states asso- 
ciated to non-linear stochastic equations exactly. 

The simplest equation, which serves as the prototype 
for our investigations, is the Kardar-Parisi-Zhang (KPZ) 
equation 



dth = vV 2 h 



(1) 



describing the non-equilibrium fluctuations of a (height) 
field h(x,t). The equation is equivalent to the Burgers 
equation (for a field u = Vh) for vorticity-free turbulence 
0| , and appears in various guises in the study of domain 
walls 0] and directed polymers 0] in a random poten- 
tial, surface growth Q, and even the gene (or protein) 
sequence alignment problem [jj ^} . Still more problems 
can be formulated as generalizations of the KPZ equation 
that accommodate spatial anisotropy or the interplay of 
several fields. Examples in the literature include the dy- 
namics of a vicinal surface [Til IT^ . the growth of two 
coupled surfaces 0] or of a magnetic surfac e IT3 . HEl Il6| . 
and the transport of a flux line or polymer [Tg • 

The stochastic aspect of Eq. Q is due to the noise 



r](x,t), which has zero mean and short-range correla- 
tions in space and time. In the absence of the non- 
linear term proportional to A, it reduces to a standard 
Langevin equation, with a Gaussian steady state. In 
one dimension, the non-linear term does not modify this 
steady state as the associated probability current in the 
Fokker-Planck equation is zero. This observation mo- 
tivates our search for other equations with this prop- 
erty, namely an easily guessed (equilibrium) steady state 
which is not affected by the additional (non-equilibrium) 
non-linearities. 

In Sec. |nl we start by constructing the Fokker-Planck 
equation for the one dimensional KPZ equation and ex- 
plicitly showing that the probability current due to the 
non-linear term does not modify the steady state, as it 
appears in an integral of a complete derivative. This ob- 
servation is then used as a basis for constructing other 
one dimensional non-linear equations that share this 
property. Indeed, we find that the class of such equations 
is quite large, including some equations already consid- 
ered in the literature. 

Higher dimensional versions of the KPZ equation may 
also obtain a Gaussian steady state in spite of their non- 
linear character. We discuss such a case in Sec. 11111 
namely an anisotropic variant of the KPZ equation in 
two dimensions with non-linear terms of opposite signs 
in the two directions. Using renormalization group meth- 
ods, Wolf [13 showed that this model indeed flows under 
renormalization to a linear fixed point. Generalizations 
of this equation with calculable steady states are also 
constructed; they all share a hidden symmetry under re- 
flection, absent in the isotropic KPZ equation. 

The examples from one and two dimensions motivate 
the search for more general principles governing the exis- 
tence of simple steady states, taken up in Sec. II VI Specif- 
ically, we consider stochastic dynamics of multiple fields 
coupled by non-linear (possibly anisotropic) generaliza- 
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tions of the KPZ term, and ask whether they admit 
Gaussian steady states. A direct solution of the Fokker- 
Planck equation becomes considerably more difficult and, 
instead, we derive two sets of general prescriptions on the 
coefficients for this to occur. These prescriptions may 
be viewed as generalized fluctuation-dissipation relations 
|18|| . and are quite restrictive. In particular, they cannot 
be satisfied in three and higher dimensions as we show in 
Sec. E| 

Having obtained general prescriptions, in Sec. I VII we 
apply them to equations for coupled fields in one and two 
dimensions. Some of the examples we discuss correspond 
to equations that have already appeared in the literature, 
in particular pertaining to the dynamics of a flux-line or 
polymer I - Sec. I VI A")l and to the growth of a magnetic film 
fSec. IVIB)l . However, in Sec. IVI CI we propose a novel 
set of equations to describe the coupling of the strain 
field of a growing crystal to its height fluctuations. Wc 
find that Gaussian steady states are indeed permitted for 
these equations in special cases. 

The Appendix treats a simple example aimed at illus- 
trating how the systematic approach may be extended to 
exactly calculable non-Gaussian steady states. 



II. GENERALIZED GROWTH EQUATIONS IN 
ONE DIMENSIONS 



Consider the probability distribution V[h] for config- 
urations of the field h(x). As the surface changes in 
time according to Eq. (Q, the corresponding probabil- 
ity evolves according to the Fokker-Planck equation 



d t V = 



v\J 2 h + -X(\Jh) 



SV S 2 V 
Sh + Sh 2 



( 2 ) 

The term in the square brackets is due to the determin- 
istic probability current and the remainder comes from 
the stochastic noise, assumed Gaussian with (??(x, t)) = 
and 



fa(x, t)T](^,t')) = 2DS d (x - x')8(t - t'). 



(3) 



In equilibrium, D = ksT; more generally, D is a measure 
of the magnitude of the noise. 

A steady-state solution is one for which d{P = 0. In 
the absence of the non-linear term, the steady-state so- 
lution of Eq. J2J is a simple Gaussian, 



V = Aexp 



2D 



d d x(V/i) 5 



(4) 



where A is a normalization constant. In general, this is 
not a steady state for A ^ 0. In one dimension, however, 
the contribution of the non-linear term to the probability 
current can be simplified to 

d t V = -V [ dx^-(d x h) 2 d xx h 



= -V 



dx- „ 
2D 

dxd- 



— (d x hf 
6D y ' 



(5) 



a surface integral safely set to zero in the limit of an 
infinite system. Thus, the steady-state spatial corre- 
lations are not influenced by the presence of the KPZ 
non-linearity and coincide with the Gaussian (Edwards- 
Wilkinson 19]) ones. We make no statements about the 
stability of the steady state of Eq. (@J ; however, here, and 
in every other example for which simulations are avail- 
able, numerical results indicate that the simple Gaussian 
steady states we discuss are indeed the ones achieved at 
long time. 

The one-dimensional KPZ equation is a particular in- 
stance of a more general class of equations, 



d t h = f (d x h) d x d x h + g (d x h) + 77, 



(6) 



with / and g arbitrary functions. These equations obtain 
a steady state 



V = Nexp 



1 

D 



dxF(d x h) 



(7) 



where / is the second derivative of F (i.e., d 2 F(u)/du 2 = 
f{u)). For g(d x h) = 0, Eq. JBJl is again a standard 
Langevin equation, while the contribution of this func- 
tion to the probability current is 

d t V = -V J dx^g(d x h)f(d x h)d x d x h 



V I dxd x 



pG(d x h) 



= 0, 



(8) 



where G is the primitive of gf (i.e., dG(u)/du — 

g{u)f(u)). ' 

The special case of a cubic non-linearity g(d x h) — 
\X(d x h) 2 + ^X'(d x h) 3 (with / = 1) was introduced to 
describe an interface separating stationary phases of the 
Toom model [2(|. The Gaussian steady state corrobo- 
rates the marginal irrelevance |2l| of the cubic term, and 
implies spatial correlations of the form 

((h(x) - Kx')) 2 ) 1 ' 2 ~ jar— x '\ 1/2 . (9) 



III. ANISOTROPIC EQUATIONS IN TWO 
DIMENSIONS 

To describe the growth of a vicinal (slightly miscut 
from a low index facet) surface, Villain introduced 
an anisotropic version of the KPZ equation, which was 
subsequently studied with a renormalization group cal- 
culation by Wolf [13]. This generalized equation has the 
form 

d t h = v x dlh + u y dlh+^\ x (d x hf + ^\ y (d y hf + r,. (10) 

Under renormalization, the subspace with v x jv y = \ x / X y 
is fixed and equivalent to the isotropic KPZ equation 
modulo a rescaling of x or y. This subspace is locally 
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attractive, so that the equation flows to a strong-coupling 
limit if X x and X y have the same sign (stability requires 
v x -,Vy > 0) . The more surprising behavior arises when 
the product A^Aj, is negative, in which case the flows 
converge to a fixed point with vanishing non-linearities. 

This vanishing of non-linearities at long length and 
time scales suggests a Gaussian steady-state probability 
density, as corroborated by an exact solution of a discrete 
model belonging to the same universality class and 
by a direct solution of the Fokker-Planck equation [23l| . 
Indeed, the Ansatz 



V = N exp 



1 

2D 



dxdy 



v x (d x h) 2 + v y (d y h) 2 



.(11) 

with the generalized fluctuation-dissipation condition 
Vxjvy — —X x /Xy, solves for the steady state. To ver- 
ify this, we note that the contributions from the non- 
linearities take the form 



d t V 



' I 2 2 

2D ' dxdy ^ x $ x ") +^v( d v h ) 



[V X d 2 x h + Vyd 2 y h] 



V_ 

2D 



dxdy 



{d x 


^X V X 


3 


+ 8y 


\yVy 

3 



(d x h) + \ y v x (d x h) (dyh) 2 



(dyh) 3 + \ x v. y (d x h) 2 (dyh) 



+2d x hdyhd x d v h (X y v x + X x v v ) } . (12) 

If \ y v x + \ x v y — 0, the above contribution is the di- 
vergence of a vector field, and hence vanishes subject to 
the usual boundary conditions. This non-perturbative 
derivation complements the renormalization group anal- 
ysis 01 which captures perturbatively the character of 
this state at large scales and the dynamics that lead to 
it. In particular, it demonstrates that a Gaussian steady 
state (and the logarithmic roughness it implies) obtains 
at any length scale, and not only in the long wavelength 
limit. 

By suitable rescalings of x and y, we can make v x = 
v y = v, so that the steady state reduces to Eq. (0} with 
d = 2. Trivially, this steady state also holds for any 
equation related to 



d t h = vV 2 h+- 



(d x h) - (dyh)' 



(13) 



by a rotation of the plane. This class comprises all the 
equations of the form 



d t h = iA7 2 /i + 



Ai r 



(d x h) 2 - (dyh) 2 



\id x hd y h + T), 
(14) 

where arctan(A2/Ai)/2 is the plane rotation angle; in par- 
ticular, the equation 



dth = vSI 2 h + Xd x hd y h + rj 



(15) 




is obtained from Eq. (|13fl by a 45° rotation. 

The surprisingly simple steady state of Eq. (|T3*)) results 
from a "hidden" symmetry under the transformation 



(16) 



The symmetry h—* —his precisely the one broken by the 
isotropic KPZ non-linearity. It is restored here, provided 
the plane is also inverted about an appropriate axis: the 
bisector in the case of Eq. i|13|) . the x- or y-axis in the 
case of Eq. <|15fl . and a properly rotated axis in the gen- 
eral case of Eq. I|14l) . This hidden symmetry sheds light 
on the renormalization group analysis |l2j . as any two 
dimensional anisotropic KPZ term may be written, upon 
rotation of the plane, as the sum of an isotropic part 
and the antisymmetric part of Eq. (|13fl whose subspace 
is invariant under renormalization. 

By analogy to Eq. ©, we can generalize Eq. (fTK|l to 
include a more complicated Laplacian term, as 

d t h = f X (d X h)d X 2 h + fy(dyh)dy 2 h + Xd X hdyh + T). (17) 

Indeed, it is easy to check that the probability density 

V = Ncxp j-i J dxdy [F x (d x h) + F y (d y h)]} , (18) 

where f x (u) = dF x (u)/du and f y (u) — dF y (u)/du, is 
stationary. 



IV. GENERAL PRESCRIPTIONS FOR 
GAUSSIAN STEADY STATES 

The solution of Fokker-Planck equations by direct 
check of Ansatze soon becomes laborious beyond simple 
one and two dimensional examples. Instead, we derive 
general prescriptions on the structure of the equations of 
evolution, for the realization of Gaussian steady states. 
We consider equations of the form 

dMx,t) = C^[h]+N^[h]+rH(x,t), (19) 

x e R d , and Gaussian 



for n coupled fields (i — 1, . . . , n 
(thermal) noise with correlator 



(» 7i (x,t>fe(x',t')> 



2 AM' 



(x-x')S(t-t'). (20) 



C x and A/" x denote linear and non-linear functionals of 
the fields, respectively, evaluated at point x. If rji and 
rjj-ti are uncorrelated [24[, the Fokker-Planck equation 
reads 



d t V 



sv_ 

Shi 



(21) 



where we have absorbed D, in a rescaling of hi (by ^/Di) , 
and reduces to 



dtV = \ d d 



-(4? 



x 'Sh t + 6h 2 



(22) 
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if C and M depend upon the derivatives of h only. We 
are looking for a Gaussian probability density 



(23) 



with Q[h] a quadratic form and N a normalizing factor, 
that solves the steady state d{P — 0, i.e. 



-(4 i) +-A4 i) ) 



SQ_ 

Shi 



SQ_ 

Shi 



(24) 

The quadratic terms cancel if C and Q are related 
through 



(25) 



and it remains to find the form of C and TV for which the 
integral 



(26) 



vanishes. This is the case if the integrand either vanishes 
identically or is the divergence of a vector field. In either 
case, the integrand is unchanged by a variation 5h(pc) 
that vanishes at infinity. If Sh is localized at x, the con- 
dition J[h + Sh] = J[h] translates, to first order in Sh, 
into 



§hj(x) 



d d yN^C^ = 0, 



(27) 



where the summation over i is understood. 

For coupled KPZ-like equations [25j of the form 



Since repeated indices are summed over, Eq. (|29|l rep- 
resents n conditions — one for each possible value of the 
index j. Each of these in fact encapsulates more than 
one constraint: Eq. (|29|l is composed of terms that come 
in one of two derivative structures and, as the equation 
must be true for arbitrary h, terms of a given derivative 
structure must cancel independently. This is ensured 
by the following two sets of generalized fluctuation- 
dissipation conditions on the tensors v and A: The 
first condition comes from grouping terms in Eq. 1|29|) 
which are products of first and third derivatives, (such 
as dadpd^hkdghi), and reads 



(!) Y,p( v ijP(a 



tjP(a)P((3)\klP( 7 )S (-ggj 
~ v ikP(a)P{fi)\jlP(i)S ) = 0, 

where the summation runs over the six permutations of 
the indices a, (3, 7. A second condition comes from group- 
ing terms of the form d a d^hkdpdshi, and gives 



( n ) J2r,R'( 2v ijR(a)R'(j)^iklR(0)R'(6) 

-VikR(a)R(J3)^ijlR' (j)R'(S) 
— VUR'(y)R'(S)\ijkR(a)R(p) 



0, 



(34) 

where the summation runs over the two permutations of 
the indices a, (3 and the two permutations of the indices 
7, 5. Each pair of conditions corresponds to a given choice 
of numerical values for j, fc, /, a, f3, 7, and S. 

Conditions (I) and (II) are necessary for a Gaussian so- 
lution of the steady-state Fokker-Planck equation. They 
are also sufficient conditions, since J\h + 5h] = Sf[h] 
to first order for any Sh implies J[h + Sh] = J[h] to 
all orders, and J^[h] = constant = c. But normalization 
of V allows only c — (otherwise V either increases or 
decreases uniformly), and consequently d{P = 0. 



d t hi = Vijapdadph + -Xijkapdahjdphk + rji, (28) 

where Latin indices i,j,k — 1, ... ,n refer to field com- 
ponents and Greek indices a, (3 — 1, . . . , d refer to spatial 
components, Eq. lf2*7|l reduces to 

Vijapkki-ysidadpdyhkdshi + d a d~fh k dpdshi) 
-vika/3^ijhs(d a d^hkdshi + d a dphkd^d s hi) = 0, 

(29) 

after some renaming of the indices and using the fact that 
v and A can always be chosen to satisfy the equalities 

Vijufi = Vijp a and XijkaP — \kjpa- (30) 
Furthermore, Eq. (|25|l requires the symmetry 

Vijafi — Vjiafi, (31) 

and the stationary probability density reads 



V = Ncxp (- J d d x'^-d a hid h. 



(32) 



ABSENCE OF GAUSSIAN STEADY STATES 
IN THREE AND HIGHER DIMENSIONS 



If the matrix of Laplacian coefficients Vij a f3 is positive 
definite, as required by infra-red stability, it can be diag- 
onalized into 



— VijS a p, 



(35) 



by successive rotations and rescalings, and the prescrip- 
tions (I) and (II) simplify correspondingly. In three 
and higher dimensions we may choose the space indices, 
in applying the prescriptions, such that a = 7 while 
a 7^ (3 S a. With this choice, it is straightforward 
to check that prescription (II) forbids a non-vanishing 
contraction of tensors with different space indices, hence 
enforces the partially diagonal form 



VsaX 



(36) 
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(where i is summed over but a is not). With this con- 
straint, prescription (II) takes the form 



J2r.,R'[ v ij^jklR{l3)R{P)8R(a)R'{-i)5R{P)R'{8) 

— ^{VikKjlR'WR'iy) + v il^ijkR(a)R(a)) 

]=0. 
(37) 

This equation expresses a set of different conditions, one 
for each choice of values of the indices that are not 
summed over. Specifically, for a particular choice in 
which a = 7 7^ j3 = S, Eq. I|37l) translates into 



^ij ^-iklc 



(38) 



For the sake of visual ease, let us define the object tp a 
VijXihiaa (the dependence of ip a on j, k, I is tacit), 
terms of which this identity reads 



(39) 



clearly, if it is possible to choose three or more distinct 
values of the indices a, /?, this condition is frustrated and 
admits only the trivial solution 







(40) 



for all j, k, I, a. Viewed as a vector identity, this requires 
that any vector Vj be orthogonal to any vector Xkiaa 
(with components labeled by i = 1, ...n). As long as the 
matrix Vij is non-degenerate, there are n non-vanishing 
independent vectors Vj (for j = 1, ...n) and Eq. I|40(l is 
satisfied only if the vectors Xkiaa vanish for all k,l,a. 
Hence, no Gaussian steady state is achievable in three 
and higher dimensions if non-linearities are present in 
the equations of evolution pr| . 



VI. COUPLED FIELDS IN ONE AND TWO 
DIMENSIONS 

In the case of a single field fluctuating in two dimen- 
sions, prescriptions (I) and (II) immediately enforce the 
form of Eq. I|14[) for which a Gaussian steady state may 
be reached, as we checked explicitely in Sec. IIIII In what 
follows, we discuss examples in which the coupling among 
several fluctuating fields broadens the class of non-linear 
equations with Gaussian steady states beyond this spe- 
cific anisotropic form (with coefficients of opposite signs). 



A. Coupled Lines and Drifting Polymers 

An array of fluctuating directed lines H, IE 01 * s 
parametrized by a single variable, consequently the Greek 
indices in Eq. I|28ll all take the same value and may be 
omitted. Equation (|28[) then reduces to 



d t hi = Vijd x d x hj + -\ijkd x hjd x hk + % (41) 



a generalization of the one dimensional KPZ equation for 
several coupled fields. In this simple case, prescriptions 
(I) and (II) are fulfilled by any tensors v, A such that 



ikMjl 



mX 



(42) 



V = N exp [ — / dx-^-d x hid x hj 



Indeed, with this Ansatz 



where Xijk can always be chosen symmetric in j, k, and 
the sum over i is understood. It is easy to check that 
these relations ensure a stationary probability density 



(43) 



d t V = - J dx-X ik id x h k d x hiVijd x d x hjV (44) 

= - J d x \ {VlMldxdx h 3dx h kdxhl 

+Vik\ijid x hjd x d x hkd x hi 
+mXikjd x h J d x h k d x d x hi)V, 

where we have renamed the mute indices to obtain the 
last equality. But Eq. (|42|l implies that all three coeffi- 
cients are identical, and 



dx\vijX ik id x (d x hjd x h k d x hi) =0. (45) 
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A special case of Eq. 

d t h\\ = v\\d x d x h\\ + \X\\{d x h\\) 2 + ^Xx(d x h±) 2 + tj\\ 
d t h± = vyd x d x h\_ + X x d x h\\d x h± + r}± 

(46) 

was introduced in Refs. 17] to describe a directed poly- 
mer drifting perpendicularly to itself. Here h\\ and h± 
are interpreted not as height fields associated with two 
different lines embedded in two dimensions, but rather as 
dynamically coupled longitudinal and transverse (to the 
average velocity of the polymer) fluctuations of a single 
line embedded in three dimensions. The conditions of 
Eq. I|42|l for a Gaussian steady state simplify to 



v\\ X 



u_ L X > 



(47) 



in agreement with a direct check 0, [2^ • 

In the simplest case with identical longitudinal and 
transverse coefficients, a stationary Gaussian distribution 
follows trivially from the steady-state properties of the 
one-dimensional KPZ equation (discussed in Sec. ITTfl . 
since the equations 



d t h\\ = vd x d x h\\ + ~X[(d x h\\) 2 + {d x hj_) 2 } +rj\\ 
d t h± = vd x d x h± + Xd x h»d x h± + r]± 



are equivalent to 



(48) 



d t h + = vd x d x h + + ±X(d x h+) 2 + r]+ , , 

d t h_ = vd x d x h_ + \X{d x h-) 2 + 77_ ' [ * J) 



G 



with h± = h\\±h±_ and r)± — n\±f]_L. Clearly, the remark 
extends to higher spatial dimensions, where 

d t h l{ = + iAKV/i,,) 2 + (V^) 2 ] + 77,1 ( 

d t h± = vV 2 h±+XVh v Vh± + ri± ' 



transform into 

d t h+ = vV 2 h + \X{Vh + ) 2 + r]+ 
d t h+ = vV 2 h_ + \\{Vh_) 2 + 7?_ 



(51) 



Therefore, the "roughness" exponent £j_ of a passive 
scalar h± advected [23,|2g by a Burgers flow u = V/i||, 
defined through 

v (^(x)-^(y)) 2 ) 1/2 ^|x-y|^, (52) 
is none other than the KPZ roughness exponent 

B. Magnetic Growth 

In a growing magnetic material, the spins may be as- 
sumed frozen in the bulk while still fluctuating on the 
surface, which itself fluctuates in height [iH llq . For 
the case of XY spins, described by a single angular field 
9(x,y,t), Ref. |16j notes that a two-dimensional version 
of Eqs. governs these non-equilibrium coupled fluc- 
tuations. In the modified notation, these equations of 
evolution read 



;n = v h V 2 h+\\ hh {Vh) 2 



d t h 



d t 9 = v g V 2 9 + X M Vh ■ V0 



+ \\gg{Ve) 2 + Vh 

+ Ve 



(53) 

(Obviously, Eqs. (|53|l fail to capture the periodic nature 
of 0, and with it the potentially relevant presence of topo- 
logical defects.) 

Using prescriptions (I) and (II), it is easy to show that 
no Gaussian steady state exists for such isotropic equa- 
tions as long as h and 8 are decoupled at the linear level. 
However, the stationary Gaussian distribution 



V = Nexp 



dxdy 



y(VM 2 + f(V^ 



(54) 



is achieved by a natural extension of the anisotropic 
Eq. O, 



d t h - v h v 2 h + ±\ hh i(d x h) 2 -(d y hy 



+ \\ee[(d x 6) 2 -(d y e) 2 } + 



Vh 



(55) 



d t e = v e y 2 e + \ h0 {d x hd x e~d y hd y e) + rie 



provided VhXgg = vgXhe (in direct analogy to Eq. I07J). 
We note also that these equations again satisfy the sym- 
metry of Eq. (fTB|> . 



C. Crystalline Growth 

The height fluctuations of a material characterized by 
internal order parameters, as in the above case of a grow- 
ing XY magnet, are subjected to the fluctuations of these 



order parameters. Conversely, the evolution of the inter- 
nal order parameters depends on the height fluctuations. 
In contrast to amorphous growth, we may say that the 
fluctuations of an ordered material results from a textu- 
ral growth, as the the additional fields invest the interface 
with a texture that constrains its fluctuations. The prime 
example is that of the growth of a crystal in which sur- 
face phonons interact with height fluctuations. In anal- 
ogy with Eqs. I|53(l . we propose the following equations 
for isotropic crystalline growth, 

" d t h = u h V 2 h + \X hh {Vh) 2 + ±aE](V • u) 2 

+ f)X ( 3l Vu, • Vui + ±Al 3 2<9iU • Vuj + i] h 
d tUi = i^W + i/^V -u + X^dihV u ' 
+\® Vh ■ Vu* + a£2 V/i • d lU + r}i 

(56) 

where u(x,y,t) is the surface displacement vector field, 

and Vu\ 
through 



and Vu \ Vu are related to the usual Lame coefficients 



(i) _ 

**u /^Lamc 



,(2) 



— /-*Lamc + A 



(57) 



Lame 



One would like to know, given the richness of Eqs. J56J, 
what phases they describe beyond the usual KPZ (amor- 
phous) phase. As a first step towards a complete answer, 
we discuss crystalline growth equations that admit an 
exact Gaussian steady state. 

Equations l|56|l tacitly assume a triangular lattice at 
the microscopic level, since other lattices would be re- 
flected in an anisotropic continuum limit that reproduces 
the appropriate crystal symmetries. In the context of 
anistropic equations, trivial extensions of Eqs. I|55|) sup- 
port a Gaussian steady state. 

Rather than dwelling on these examples, we turn to a 
new possibility, namely isotropic equations of the form 
of Eqs. i|56|l that admit a stationary Gaussian probabil- 
ity density, as allowed by the presence of linear couplings 
between u x and u y . A somewhat tedious but straightfor- 
ward application of prescriptions (I) and (II) to Eqs. I|56|) 
yields two (and only two) non-trivial solutions, charac- 
terized by a vanishing shear modulus iffi — and a 



(i) 



,(2) 



vanishing bulk modulus i>u 

With a vanishing shear modulus v^ 1 = 0, (I) and (II) 
imply 



0. 

,(!) 



A;,/, — a£5 — aI5 = a! 2 ? — a1 3 ' 







and 



^a« = 4 2) a« 

resulting in the equations of motion 



d t h = v h V 2 h+ §A„ U (V -u) 2 + i] h 
d t u = u u V(V ■ u) + X hu Vh(V ■ u) + r) 



(58) 



(59) 



(60) 
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where the coefficients have been renamed (v u \ 
Ku = A-Lu, X hu = ^hu)- Tne probability density 



V = iVexp 



dxdy 



y(V^) 2 + ^(V.u) 



is stationary, as the direct check 

1 



d t V 



dxdy 



-X UU (V ■ n)v h V z h 



A uu V/i(V ■ u)i/ u V(V • u) 



dxdy 



-V ■ [Vh(V ■ u) 2 ] 



(61) 
(62) 

(63) 
(64) 

(65) 

(66) 

we may interpret Eqs. (|6U[) as describing a liquid that is 
flowing on a fluctuating surface h(x, y, t) via the partially 
decoupled set of equations 

d t h = v h V 2 h + \X uu p 2 + r) h 
dtp = v u V 2 p + A ft „V • (pVh) + V ■ t] (67) 
d t n = d x 7] y - d v r\ x 

This is again reminiscent of the advection of a scalar 
(density fluctuations are advected along height 
gradients) which is not quite passive, as p influences the 
evolution of h. Equations l|67|) describe also a special case 
of the coupled growth of a binary film or, equivalently, of 
a scalar (Ising) magnet |l5j . 



= 0, 

confirms. In terms of the density fluctuations 
p = V • u 

and the vorticity 



With a vanishing bulk modulus v. 
(II) imply 



(i) 



.(2) 



0, (I) and 



Xhh 

and 



AW 

UU 



X 



(1) 

/),-(./. 



A (2) 

UU 



A (2) 



AS=A^+A^=0 (68) 



U u A hu ' 



(69) 



resulting in the equations of motion 



d t h 



where v u 



vhS7 2 h + \X uu diUj (diUj - djUi) + r\ h 
u u dj (djiit - d t Uj) , 
+X hu djh (djUi ~ diUj) + T) 

(70) 

Introducing 



, A 



= \W 



mi = X U u , and Xhu — A^, (j 



a fictitious direction z according to the natural definition 
V x u = (d x u y — d y u x )z with V = (d x ,d y ,0) and u = 
(u x , u y ,0), we can rewrite Eqs. (|70|l in the more compact 
form 



d t h = v h V 2 h + |A„ U (V x u) 2 
d t u = —v u V x (V x u) 

— A/mVft x (V x u) + rj 



Vh 



(71) 



With this form at hand, the stationarity of the probabil- 
ity density 



iVexp 



dxdy 



V}, 



(VA) a + £(Vxu) s 



. (72) 

is obtained in direct analogy with the case of vanishing 
shear modulus. In terms of p and f2, Eqs. I|71|l become 



d t h = v h \l 2 h 
d t il = v u V 2 Vt + X 
dtp = V\\ 



hu v ■ (nv&) 



(73) 



none other than Eqs. I|67l) with p and Q interchanged, 
but still a transverse noise r/± = d x r\ y — d y r/ x driving the 
fluctuations of ft and a longitudinal noise 7711 = d x r\ x + 
d y r\y driving p. Here we can interpret Q as the density 
of dislocations whose presence locally affects growth. 

Experiments on surface growth do not point to a 
unique characterization |30| and, in particular, most mea- 
sured roughness exponents differ from those expected on 
the basis of the KPZ equation |sL l3ll . This may re- 
sult from conservation laws, incorporated in some molec- 
ular beam epitaxy models |33l IsiL I35L l36T |. but it also 
may be the consequence of the dynamic coupling of the 
height fluctuations with the fluctuations of the intrinsic 
order parameter of the material, as in magnetic or crys- 
talline growth. Thus, having established that the special 
Eqs. (|60|) and (|71|) admit Gaussian steady states, one 
would like to know whether the full Eqs. I|56l) may flow to 
them under renormalization. If Eqs. I|6U|I or (|71|) indeed 
have a basin of attraction, we are left with the surprising 
conclusion that the coupling to crystal vibration teth- 
ers the fluctuating surface, in an appropriately prepared 
sample, to logarithmic roughness, in marked contrast to 
the KPZ roughness associated with amorphous growth. 

Rather than renormalizing the full set of Eqs. (|56|l . one 
may, as a first attempt, consider the renormalization of 
Eqs. (E3 or lO with an added KPZ term {X hh ^ 0), 



v = v h V 2 h+ \x hh iyh) 2 + \X uuP 2 



d t h 

d t p = v u V 2 p + X hu V ■ (pVh) + r?|| 
<9 t Sl = rjA 



Vh 



(74) 



The sub-space of Eqs. I|74(l is closed under renormaliza- 
tion since they are invariant under the transformation 

u^u + VxA, (75) 

with A(a;, y) an arbitrary smooth function. (Similarly, 
Eqs. (|73|l are invariant under 



V0, 



(76) 



with <f> arbitrary.) Also, Eqs. I|74|l are interesting in their 
own right: the field p may be interpreted as density fluc- 
tuations, e.g. of surfactants, sliding on the surface. Sim- 
ilarly, if h describes the height of a material surface (such 
as a liquid film) , p may play the role of material surface 
density fluctuations (proportional to the thickness of the 
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film in the case of an incompressible liquid). These den- 
sity fluctuations are conserved as the particles, or liquid 
elements, slide along gradients in the surface height, ac- 
cording to Eqs. I|74l) . However, the implicit p — > — p 
symmetry of Eqs. (|74|l is difficult to justify physically 
for surfactants or liquid films. 

As a final non-perturbative observation on Eq. (|74|l . 
we note that they are invariant under the infinitesimal 
tilt 



x + Xet 
h + e ■ x 



(77) 



in the sub-space Xhh — Xhu = A, so that Xhh and Xh u are 
not modified by coarse-graining. Their renormalization 
flows are then given by 



(78) 



where Zh and z p are the dynamical exponents associated 
to the fields h and p, respectively, and (h is the height 
roughness exponent. If a fixed point occurs at Xhh = 
A/m 7^ 0, the relation 



Zh + Ch = z p + Ch = 2 



(79) 



holds exactly, while stability of a fixed point with Xhh = 
Xhu = requires 



Zh + Ch < 2 and z p + Ch < 2. 
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APPENDIX A: NON-GAUSSIAN STEADY 
STATES 

In this Appendix, we show through a simple example 
how the method used in the bulk of the paper may be 
extended to search for non-Gaussian steady states. We 
focus on equations of the form of Eq. (|17fl in Sec. II I II 
where the usual Laplacian smoothing term is promoted 
to a more general object. 

In parallel, we promote Q (see Eq. (J221) to include 
higher order terms, as 



Q 



d d x (vijapdahtdphj 



+TTijkaf3yd a hidphjd~ ( hk 

+Pijkiai3~f5d a hid0h j d-fh k dshi H ) , (Al) 

while still imposing the relation of Eq. I|25(l between C 
and Q. Following the procedure of Sec. llVl in the special 
case of a single field and in two dimensions (for the sake 
of simplicity), we obtain the additional prescriptions 

71 ' xxQ^Xyy ^^xyaXxy -\- 7Tyy a X X x 0; (-^2) 

for any a, and 

Pxxafi^yy ^PxyaftXxy ~\~ Pyya.^X X x — 0; v^"^) 

for any a, (3. Here it and p have been symmetrized, and 
a, (3 denote x or y. 

Let us focus on the quartic term with coefficient p. For 
rotationally symmetric tensors 



and 



X a (3 = X5 a (j 



(A4) 



(A5) 



Eq. (| A3|) has only the trivial solution A = 0. By contrast, 
an antisymmetric non-linear term with 



X X x — ~^yy — A, and X xy — 0, 
allows the non-trivial solution 



Pxxxx — Pxxyy — Pyyyy — Pi 
Pxxxy — Pyyyx — P2/2 



(A6) 



(A7) 



with arbitrary p\ and p2- Thus, the equation of evolution 

d t h = vV 2 h + pi[(Vh) 2 V 2 h + d x hdyhd x d y h] 
+p 2 [d x hd v hV 2 h + [Vhfdxdyh] 

+ ^[(dxh) 2 -(dyh) 2 ]+ v (A8) 
obtains the non-Gaussian steady-state distribution 

Quartic = iVexp^- J dxdyfyvhf 

+ ^[(dxh) 4 + Q{d x h) 2 {d y h) 2 + [dyhf] 
+ j[(d x h) 3 d y h + {d v hfd x h])) . (A9) 

We note that this equation of evolution, too, satisfies the 
"hidden" symmetry of Eq. l(TH|l . 
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